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Abstract 

Interacting dark energy and the holographic principle offer a possible way of addressing the 
cosmic coincidence problem as well as accounting for the size of the dark energy component. The 
equilibrium points of the Friedmann equations which govern the evolution behavior of dark energy, 
matter, and curvature components can determine the qualitative behavior of the cosmological 
models. These possible equilibria and their behavior are examined in a general framework, and 
some illustrative examples are presented. 
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I. INTRODUCTION 



The study of type la supernovae indicates that we hve in an accelerating Universe. This 
acceleration can be described by adding another component to the cosmic inventory that 
behaves in a similar way to the famous cosmological constant. Observations are consistent 
with about 70% of the density of the nearly flat universe is in the form of dark energy with 
the remaining 30% of ingredients is composed of matter, baryonic or non-baryonic. During 
the expansion of the Universe, one component tends to replace another in a small number 
of e-foldings, so the current closeness of these two densities is remarkable. This is known as 
the cos,„,c coi,.ide.ce prob.e. Q. A perhaps assoe.ated quest.on — why the da* 
energy density is so small compared to the other physical scales in particle physics These 
issues have led to the consideration of models that are more general than adding a simple 
cosmological constant. One hope is that a model based on a simple, single physical principle 
can account for both the coincidence and the observed size of the dark energy density. 

A class of models that has been employed to explain the coincidence problem involves a 
nonzero interaction between dark energy and matter. In these theories, dark energy converts 
into matter in such a way as to create an equilibrium balance at late times. The choice of an 
appropriate interaction can effect a smooth and gradual transition from a matter-dominated 
(decelerating) Universe to a dark energy dominated (accelerating) one. These models lack a 
fundamental explanation for the source of the interaction. They are nevertheless an attractive 
way to investigate, at least at the phenomenological level, ways to account for the observed 
cosmological quantities such as the dark energy density parameter Q\ = 0.7 and its equation 
of state wa ~ — 1. If dark energy decays into matter, it is possible to imagine an equilibrium 
situation in which the tendency of the dark energy component to dominate as the Universe 
expands, is balanced by its decay into matter. These models represent deviations from, and 
are more general than, the model with a cosmological constant (ACDM). 

The problem of the scale (size) of the dark energy component has led to the consideration 
of the holographic principle. This idea offers a possible mechanism for tying the scale of 
the dark energy to the geometric average of the Planck scale and the size of the Universe. 
The authors previously have used interacting dark energy in the scheme of holographic 
principle to find a relation between the equations governing the dark energy and matter 
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densities in a flat universe \3]- A holographic principle relates the density of dark energy 
to some horizon. Some examples are the apparent Hubble, future event horizon or particle 
horizon [ill, Q, Q. It has been previously noted that, without any 

interaction between dark energy and matter, only the future event horizon gives results 
in agreement with observations 13]. However, if we introduce an interaction, all of the 
above choices for the horizon might be consistent. There exists, for example, a model in 
which a constant interaction with Hubble horizon in a flat universe gives rise to the desired 

n 

The consideration of nonzero curvature is the next obvious generalization [16| . Although 
we know from the cosmic microwave background that we live in an almost flat universe, this 
may not have been the case at all times. Furthermore, the decay of dark energy to matter 
involves the conversion of a component which tends to make the universe more flat with one 
that has the opposite behavior. Nonzero curvature also makes the evolution equations more 
involved, and one can accommodate additional flxed point solutions. 

We describe the equations that govern the dynamics of dark energy and curvature in 
Section II. Assuming that the radiation component is negligible, these are the only two 
independent density parameters out of the three (matter component, dark energy component, 
and the "curvature component") related by Qm + = I + ^k- To solve the evolution 
equations one needs to know the effective equations of state for dark energy and matter. 
This requires a speciflcation of the holographic condition on dark energy, i.e. whether the 
length scale associated with the dark energy density is the Hubble horizon, the particle 
horizon, the future horizon, or some other physical scale. In addition, one needs to specify 
some interaction. With these two speciflcations, solving the resulting differential equations 
analytically is in general impossible, and solutions have to be done numerically. Therefore, 
it is desirable to arrive at some general characterization of the types of asymptotic solutions 
that can be obtained. One needs to flnd the steady-state solutions of the equations, where 
the variable parameters are flxed for large times in the future or past. These equilibrium 
solutions, and the conditions giving rise to them, are discussed in Section HI. Some of these 
conditions have been exploited previously in the literature. To further clarify the main 
points of our results, we describe different examples in detail in the Section IV. It is worth 
mentioning that interacting models are not limited to models which employ an interaction 
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between dark energy and matter. There exist models in which there are two components 
for dark energy: a variable cosmological constant, A, and a dynamical "Cosmon", X, with 



possible interaction between them. These models, called AXCDA/ 



of addressing the cosmological constant problem. See e.g. 



13, 



models, are also capable 



19|. 



II. DYNAMICAL EQUATIONS FOR AND flk 

We start by assuming that the universe consists of matter and dark energy, nonzero 
curvature, and with negligible radiation. The evolution equations for dark energy and matter 
are, 

Pa + 3H{1 + wa)pa = -Q , 

We shall refer to the quantities Wm and wa as the native equations of state. The appearance 
of the same Q in the two equations insures the overall conservation of the energy-momentum 
tensor. Positive Q can be interpreted as a transfer from the dark energy component to the 
matter component. Presumably this interaction could arise from some microscopic mech- 
anism and might be specified more fully when the nature of the dark energy is known. It 
can be noted, however, that the consideration of quantum fields in curved backgrounds gives 
rise to particle production. However, in a cosmological scenario based on the Friedmann- 
Lemaitre-Robertson- Walker (FLRW) metric, the production rate is only at the level of the 
Hubble constant (and not the geometric average of the Hubble scale and the Planck scale). 
It is interesting to note that the required physical scale of the interaction involves, like the 
dark energy itself, the size of the universe. This suggests a holographic interpretation for 
any such interaction. 

The quantity pm is the matter component for which one usually takes Wm = 0. The 
quantity pA, the density of dark energy can be related to some length scale via a holographic 
principle. We define this length scale La according to the equation 

Here c is a constant of order one (inserted here to preserve agreement with the literature) 
and Mpi is the Planck mass. Various choices considered for this length scale are the Hubble 
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horizon Rh = 1/H, the particle horizon (PH) defined by 

t 



RpH — a 



and the future event horizon (FH) defined by 

oo 

f dt 



da 



da 



(3) 



(4) 



where a{t) is the scale factor in the FLRW metric. 

The interaction can be combined with the native equations of state, wa and Wm, to define 
effective equations of state as 



PA + 3H{1 + wf)pA = , 

Pm + SH{1 + W'^)pm = . 



(5) 



Defining the ratio r = Pm/pA and the rate T = Q/ Pa and assuming a pressureless matter 
component, for which Wm — 0, these effective equations of states are 



w 



cff 



w, 



cff 



The density parameters 



Sttpa 



Snpn 



1 r 

r3^ 



H^a^ 



are defined so that the Friedmann equation, 

SttG 



k 



(PA + Pm) - ^ ■ 



(6) 



(7) 



(8) 



gives 



The ratio r is related to density parameters by 

1 — Qa + ilk 



r — 



and its time evolution is 



r = 2>Hr 



WA-Wm + 



1 + r r 



r 3H 



3Hr 



wf - 



(9) 



(10) 



(11) 



It is evident from Eq. flTT]) that the condition w'^ = w'^ involving the effective equations of 
state gives r = 0. When this equihbrium occurs, the ratio of dark energy and dark matter 
densities is a constant. The Friedmann equation, Eq. ([8]), and the time evolution of densities 
can be combined to find the time evolution of Hubble parameter, 

IdH 3 1^ 3 ^ 

"o ~ o^'fc ~ TTWAiiA , (12) 



H dx 2 2 2 

where x = ln{a/aQ) with some fixed scale factor a^. The density parameters satisfy the 
differential equations, 

dVLp^ 



dx 
dflk 
dx 



1 r 



n^ii + + 3wAnA) . (is) 



It was shown in Ref. [3| that it is sufficient to have two physical conditions to determine the 
evolution of the parameters of the universe. The possibilities include a holographic condition 
on the dark energy p^, an assumption about the nature of dark energy and its native equation 
of state, w\ or an assumption for the form of interaction Q, or equivalently F. These three 
specifications are related by 

T = 3H{-1-wa) + 2^. (14) 

This equation demonstrates that the interaction should generically be of the same size of the 
Hubble parameter, and also suggests that holographic definitions for the interaction may be 
useful. 

Since the physical interpretation of the effective equation of state is clear, one can elimi- 
nate wa and F in favor of wj^ and to obtain the equations, 

dQjf 



dx 

"^^^ - 3^]fc^]A« - <^) + ^]fe(l + + 3<^) . (is) 



dx 
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and with the substitution 
2|. The location of the fixed 



These equations are consistent with the analysis of Ref. 
f2fc = 0, one recovers the equation for the flat case from Ref. 
points and equihbria of these coupled differential equations will determine their asymptotic 
behavior. The appearance of the factors w^^ — w'^ and 1 + are easy to understand on 
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physical grounds. The first factor merely compares whether dark energy or matter comes 
to dominate as the universe expands. The second factor compares matter to curvature 
("w^^ = — I")) so it measures whether the density of matter increases or decreases as the 
universe expands. 



III. POSSIBLE EQUILIBRIUM SOLUTIONS OF THE EVOLUTION EQUATIONS 

Equilibrium solutions occur when the right hand sides of the differential equations, 
Eqs. f fTSj) . vanish. We wish to distinguish different kinds of equilibria. First, there are 
fixed points in Q^. space. These are universal in the sense that they will be present 
unless the effective equations of state have singularities as a density parameter vanishes. 
Second, some equilibria consist of one constraint on the effective equations of state and a 
fixed value of Qa or Qk- This represents a fixed point which depends in detail on the specific 
form of the interaction and holographic definition of the dark energy. Finally, there is an 
equilibrium solution which is governed by just constraints on effective equation of states. 

A procedure to identify these equilibria is the following: Set the derivatives in Eqs. f[T^ 
to zero. Taking linear combinations of the right-hand-sides, one finds solutions 



f]A«-0 = 0, 

+ 3w'^) = . (16) 

One obtains solutions by identifying points where the density parameters vanish and/or 
where the effective equations of state satisfy a certain condition, namely, wj^ = or 
1 + 3wf = 0. 

The behavior near these fixed points is important in understanding the qualitative be- 
havior of the evolution. The fixed points can be categorized as either repellers, attractors or 
saddle points. Consider the two differential equations in the following form 20|, l2l|, |22| , 

^ = m) ■ (17) 

ax 

In this equation ft = Qk) and f (O) represents the right hand sides of the two differential 
equations of the evolutions of the density parameters. The density parameter for matter, 
flm is then given by Eq. (Q. The behavior at a specific equilibrium point depends on the 
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eigenvalues of the 2-by-2 matrix A = If the eigenvalues of this matrix are all positive, 
then we have a repeller fixed point. In other words, as the Universe expands and x increases, 
the Universe will tend to move away from the fixed point. If the eigenvalues are all negative, 
we have an attractor fixed point. Finally if one eigenvalue is positive and one is negative, 
the fixed point behaves hke a saddle point. 

There are three fixed points for the system of equations which occur at certain values 
of the density parameters. The character of these fixed points is easy to identify from the 
properties of the effective equations of state. The first fixed point occurs at f2 = (0,-1). 
This is a negatively curved solution, and would be a future asymptotic state of the Universe 
if both > — I and wj^ > — |. Since the Universe appears flat, we know that it is unlikely 
to be a viable solution. The second fixed point occurs at f2 = (0,0). This corresponds to 
a universe filled with only matter. This could be the final state of the Universe, provided 
wf^ < — I (so that the matter density grows as the Universe expands) and that wj^ > 
(so that matter density comes to dominate over the dark energy density). Finally, the de 
Sitter Universe is described by O = (1,0), a universe just filled with dark energy. This 
solutions pertains to the future asymptotic state of the Universe provided wj^ > (so that 
the matter density grows as the universe expands) and w'j^ < w'^. The eigenvalues of the 
matrix A are shown in Table [B If one replaces the effective equations of state with their 
native values {wm = for matter and W\ = —1 for a cosmological constant), one recovers 
the usual evolution behavior associated with the ACDM. One recognizes from Table [T] that 
the behavior in a more general context is governed entirely by the values of the effective 
equations of state at the fixed point. 

It is also worth mentioning that, for specific interactions, some of these fixed points can 
be eliminated. For example, if the interaction produces a w'^ proportional to something 
sufficiently singular like I/^a, then the point ft = (0, 0) will cease to be a fixed point as can 
be seen in Eq. (Q. 

The behavior of the evolution near these fixed points depends on the effective equations 
of state in Table [H the different possibilities for their dynamical behaviors are shown in 
Table [TTl [ml and IVl 

The fixed points just discussed are the familiar ones that result from behavior of each 
density component as the universe expands. There are also fixed point, or equilibrium 



8 



fixed point 


eigenvalues of A 


f2 = (0,-l) 


-(l + 3t(;f ) 

-(l + 3<^) 


n = (0,0) 


-3« - wf) 
l + 3wf 


O = (l,0) 

nm = o 


3{wf-wf) 
1 + 3wf 



TABLE I: Equilibria of the universal fixed points in fl space. The signs of the eigenvalues are 
given by comparing the effective equations of state to the that of the curvature component (which 
is —1/3), or by comparing them to each other. 





^ 3 


> 3 


y^eff ^ _1 
"^m ^ 3 


repeller 


saddle point 


^eS ^ _1 


saddle point 


attractor 



TABLE IL Characteristic behavior of the fixed point = (0,-1), when = 0. 

solutions, that can result because the effective equations of state satisfy some constraint. 
There are two equilibrium solutions where just one of the density parameters is determined 

and the other condition is a constraint on the equations of state. The first occurs when 
= and wf = — |. This equilibrium would be possible only if the chosen interaction 
gives rise to ^ ~ tt- around = 0. Although this doesn't directly make the interaction Q 





wf < -i 


wf > -I 


'^m ^ 3 


attractor if wf > wf 
saddle point if wf < wf 


attractor 




repeller 


repeller if wf < wf 
saddle point if wf > wf 



TABLE IIL Characteristic behavior of fixed point = (0,0), when 0^ = 1. 
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„,,cff ^ 1 
< 3 


„,,cff ^ 1 
'^A > 3 


^neff < -1 
"^m ^ 3 


saddle point if wj^ > wf^ 
attractor if w'^ > w^^^ 


repeller 


"'m ^ 3 


attractor 


saddle point if wj^ < w'^ 
repeller if wj^ > 



TABLE IV: Characteristic b ehavior of fixed point ft = (1,0), when 0^ = 0. 



fixed density parameter 


constraint on effective EOS 


eigenvalues of A 


17a = 


y^eff _ _1 
"^m 3 


-(1 + 30 


^k = 


eff eff 
Wm = W^A 


-317A(l-f7A)(5^-^) 
1 + 3w;f 



TABLE V: Equilibria specified by one fixed density parameter and one constraint. 



infinite, it forces F to be infinite around Q/y = 0, which might be problematic (see Eq. ([6])). 
The second equilibrium occurs when fi^ = and = wj^. This solution has been 



employed [3|, [lO|, |23j to represent the final asymptotic state of the universe, and it has to 
be an attractor to be a suitable candidate. This implies that f = from Eq. (fTTj) and the 
density parameters remain constant, provided this is a stable equilibrium as determined by 
the eigenvalues at the fixed point. The current epoch of the universe may be nearing such a 
fixed point. On the other hand, it is also possible that the current state of the universe is 
a transient one in the sense that behavior near the equilibrium point behaves like a saddle 
point. The eigenvalues of A are listed in Table IVl 

The last equilibrium condition occurs when the constraint = = — | holds. The 
equilibrium solution is easy to understand as both effective equations of state scale like 
curvature, so that the dark energy density and matter density parameters are constant. 
It is harder to see the behavior of this equilibrium as the derivative matrix. A, includes 
complicated function of the derivatives of the effective equations of state. An example of a 



solution satisfying these constraints was presented in Ref. 



15|. There the model consisted 
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of a holographic length scale set equal to the Hubble parameter, and a constant interaction 

r. 

These examples indicate the large variety of behavior that can result with fixed points. 
The behavior of the equations at these other equilibria is more involved than the universal 
fixed points in Table [H as the eigenvalues of A depend on the derivatives of the effective 
equations of state as well. Some examples will be considered in next section to illustrate the 
different possibilities. 

IV. EXAMPLES 

In this section the behavior of the density parameters will be illustrated with some choices 
for the interaction. As mentioned earlier, two conditions are required to completely specify 
the equations. If one of the conditions is determined by a holographic principle, it will 
restrict the possible behaviors of the equilibria. If one chooses the length scale entering the 
equation of density of dark energy to be the future event horizon, defined by Eq. (jll), then 
it can be shown that the effective equation of state for dark energy can be written as 24 1 



On the other hand, if the length scale is chosen to be the particle horizon, defined by Eq. (j3]), 
with a similar procedure one can show that 



In general in the FH case, wj^ < -i, and for the PH case, > —3. Therefore, after 
specification of this condition, the behavior of the fixed points will be determined by the 
behavior of , which is in turn determined by the form of the interaction. 

For all the examples in this section, the holographic length scale, entering in the equation 
of dark energy density, will be taken to be the future event horizon. La = Rfh defined in 
Eq. (jll). The effective equation of state for dark energy is given by Eq. ( ITSl) . What remains 
is a specification of the form of the interaction. The simplest choice is no interaction at 
all, for which w'^ = 0. In this case, there exist just the three equilibrium points listed 
in Table [H The behavior of the fixed points can be obtained by looking at Tables [Tll IHH 




(18) 




(19) 
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Qa 




FIG. 1: (a) Density parameters in a model with no interaction, (b) The associated flow diagram 
in the J7a, plane. The solid line represents the solution for the initial condition assumed in (a). 

and IIVI considering the fact that = and wf < -|. The fixed points at ft = (0^, 
ft = (0, —1), and O = (1, 0) are a repeller, a saddle point, and an attractor respectively [2 61]. 
The behavior of the density parameters are shown in Fig. [1], where an initial state is chosen 
so that the Universe is filled with matter, while dark energy and curvature are very small. 
The interaction drives the Universe toward a final equilibrium in which the Universe is filled 
with dark energy. The Universe experiences a transition toward and then away from a saddle 
point where it is negatively curved. The physical interpretation of this behavior is similar to 
but not the same as ACDM: matter tends to create curvature in the universe, but eventually 
the dark energy component comes to dominate and reduces the curvature. 



The next example employs an interaction between dark energy and dark matter and 
illustrates the appearance of a new equilibrium point. In terms of F, the interaction is 

r b'^n + Qk) 



3H Qa"" 

which has been used in some recent works to explain the cosmic coincidence problem 



(20) 
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25]. 



When n = 1, the interaction term Q = p/J? is just dependent on H, the Hubble parameter. 
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FIG. 2: (a) Density parameter evolution for the interaction in Eq. ()20p . when 6^ = 0.26. (b) 
Flow diagram for the evoloution of density parameters using interaction Eq. (j20|) . The solid line 
represents the solution depicted in (a), (c) Effective equations of state for dark energy and dark 
matter. 

The density parameters, acting under this interaction, approach an equihbrium solution 
similar to their observed values at the present time. The effective equation of state for 
matter is 

= -^^(1 + ^^)^a'" . (21) 

1 -|- — 

It can be seen from Eqs. f|T5l) and the definition of w'^ that this interaction eliminates 
ft = (1,0) as an equilibrium. However, with an appropriate choice for 6^, it produces a 
new attractor equilibrium at f2 = (0.7,0). This new fixed point arises from an equilibrium 
specified by = w'f as in Table El This attractor can be employed to represent the 
future asymptotic state of the Universe. Although there is no theoretical justification for 
this interaction, it turns out to be very useful in understanding how introducing interactions 
can change the behavior of the density parameters. The initial values for density parameters 
don't have any effect on their asymptotic behavior. Initial values, though, dictate how the 
density parameters evolve before reaching their final asymptotic value and how near they 
approach the saddle point-hke equilibria. 



In Fig. [2], the Universe, is taken to be filled with matter and to have very small curvature 
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initially. It moves toward a saddle point where the curvature is almost —1 and then evolves 
toward its final situation, where ~ 0.7. How close the Universe approaches the interme- 
diate saddle point depends on the initial values for Q\ and Qk- If those initial conditions are 
chosen to be smaller than the ones chosen in Fig. [2], the evolution turns around at smaller 
absolute values of curvature density parameter and then proceeds toward its asymptotic 
solution. For the last fixed point, the effective equations of states are equal, and the matter 
and dark energy density parameters remain at their equilibrium solutions. Any curvature is 
eliminated as the Universe expands under the infiuence of the presence of the dark energy 
component. This generalizes our previous result where an exactly fiat initial condition was 
assumed j^. 

A final example, involves an interaction which causes the universe to explore a transient 
state for a long period of time by exploiting a fixed point which is a saddle point. Weighting 
the interaction, Eq. fl2U]) . by a factor of exp(— pfi^), we have 

3H f^A" ■ 

The exponential factor, while somewhat artificial, guarantees that evolution is driven close 
to a saddle point. The effective equation of state for matter is 



eff -fe^e-^'^Hl + ^fe)^A 



n 



We do not ascribe any physical motivation for such an interaction, but merely employ it 
to show how the evolution can be driven close to a saddle point. We again choose n = 1 
for the numerical solutions. In this case, the Universe starts from a repeller, approaches 
the saddle point, then is finally repelled toward the attractor. In Fig. [3] the initial value of 
the curvature is set very close to zero and negative, while in Fig. HJ the initial value of the 
curvature density parameter is set very close to —1. In both cases, however, the Universe is 
pushed by the interaction toward the saddle point where we have fi^ ~ —0.05 {27], then goes 
toward the attractor located at fi = (0.7,0). At this point, again the condition = w*]^ 
holds, as can be seen in Fig. H] (c). 
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FIG. 3: (a) Density parameters in an interacting model when F is proportional to 6"^^*=, when 
6^ = 0.26 and p = 5. (b) The initial value of the curvature density parameter is taken to be a very 
small and negative number. 




FIG. 4: (a) Density parameters in an interacting model with interaction, Eq. (j22p . The initial value 
of is set close to —1. The interaction drives toward zero very rapidly, (b) The associated 
flow diagram in the J^Aj^^/c plane. The solid line represents the solution depicted in (a), for the 
region Qk > —0.1. (c) Effective equations of state for dark energy and dark matter. 
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V. CONCLUSION 



The behavior of the Universe, where there is an interaction between dark energy and 
dark matter, has been studied in the presence of curvature. In the hterature, many different 
scenarios, obtained by choosing different interactions and different choices for length scale 
which enters the definition of dark energy density, have been explored. The present work 
unifies these attempts in this approach toward solving cosmic coinsidence problem, at the 
qualitative level, by studying the possible equilibria showing up in the dynamical equations 
of dark energy and dark matter. 

Understanding the fixed point equilibria is important, because the asymptotic behavior 
at these points is insensitive to initial conditions. Attempts to solve the coincidence problem 
in the hterature usually focus on stable solutions for which Q). — and — w^. With 
appropriate choice of interaction, a fixed point can occur at ft\ ~ 0.7 and Jl^ ~ 0.3, the 
observed values for these parameters. We showed that adding curvature to cosmic inventory 
results in additional equilibium solutions. For example, the equilibrium point Q = (0,-1), 
which corresponds to flrn = 0, doesn't exist when there is no curvature. The new and old 
equilibrium solutions, and the possible behavior of them, are studied and categorized in 
Section III. A particular theory with an interaction between dark energy and dark matter 
doesn't typically possess all of these equilibria. For example if there is no interaction at all, 
there exists just the three equilibria listed in Table I. By a suitable choice of interaction, 
one is able to eliminate some of the equilibria and add new ones. The examples provided in 
Section IV illustrate how this process works. 

Observation indicates that there is little curvature at the present epoch. Additional 
equilibria can be obtained by considering the Universe at all eras, where curvature is not 
necessarily small. Any saddle point equilibrium could represent a transient era (but perhaps 
long-lived) of Universe. Any attractor equilibrium, providing initial conditions let Universe 
be driven toward it, would represent the final state. In fact, no observation has confirmed 
that Universe is in a stable solution of its dynamical equations. So it is quite possible that 
Universe is currently experiencing a saddle point-like equilibrium and will ultimately be 
deriven toward a fixed point with complete different properties. 

The supernovae data requires a recent transition from deceleration to acceleration, so a 
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consistent solution requires that we are arriving at this equihbrium at the present epoch. 
The fixed point solutions represent an amelioration of the rapid transition to dark energy 
domination such as occurs in the ACDM. Any particular choice of interaction could be 
supported or ruled out by comparing its results with recorded history of Universe. The 
examples provide insight into the classification of the equilibrium points. 
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